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Abstract: We extend the analysis of of the Standard Model Higgs inflation accounting 
for two-loop radiative corrections to the effective potential. As was expected, higher loop 
effects result in some modification of the interval for allowed Higgs masses m^^^ < rriH < 
^max, which somewhat exceeds the region in which the Standard Model can be considered 
as a viable effective field theory all the way up to the Planck scale. The dependence 
of the index of scalar perturbations on the Higgs mass is computed in two different 
renormalization procedures, associated with the Einstein (I) and Jordan (H) frames. In 
the procedure I the predictions of the spectral index of scalar fluctuations and of the tensor- 
to-scalar ratio practically do not depend on the Higgs mass within the admitted region and 
are equal to Ug = 0.97 and r = 0.0034 respectively. In the procedure II the index acquires 
the visible dependence on the Higgs mass and and goes out of the admitted interval at mn 
below mmin- We compare our findings with the results of M. 
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1. Introduction 

A theory of inflation provides an attractive explanation of the basic properties of the Uni- 
verse, including its flatness, homogeneity and isotropy (for a recent review see [^). In 
addition, it gives a natural source for generation of an almost flat spectrum of fluctuations, 
leading to the structure formation (for a comprehensive discussion see [Q]). Many realisa- 
tions of inflationary scenario introduce new physics between the electroweak and Planck 
scales, and postulate the existence of a special scalar field — inflaton. Different models 
of inflation lead to distinct predictions of the main inflationary parameters, such as the 
spectral index Ug of the scalar perturbations and the tensor-to-scalar ratio r. A number of 
models have been already excluded or are in tension with the cosmological observations of 
the cosmic microwave background (CMB) 

Within the variety of inflationary models there is one which plays a special role. It 
does not require introduction of any new physics and identifies the inflaton with the Higgs 
field of the Standard Model (SM) Q. The key observation which allows such a relation is 
associated with a possible non-minimal coupling of the Higgs field H to the gravitational 
Ricci scalar R, 

-^non— minimal H R . (f 
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For large Higgs backgrounds ^/i^ > Mp (here Mp = 2.4 x 10^^ GeV is the Planck scale and 
= 2H^H) the masses of all the SM particles and the induced Planck mass [Mj^^]^ = 
Mp + ^/i^ are proportional to one and the same parameter, leading to independence of 
physical effects on the magnitude of h. In other words, the Higgs potential in the large- 
field region is effectively flat and can result in successful inflation. The fact that the non- 
minimal coupling of the inflaton to Ricci scalar relaxes the requirement for the smallness 
of the quartic scalar self- interaction was noted already in 0, 10, 11, 12, p^ ]. 

This model of inflation is very conservative in a sense that it does not require the ex- 
istence of new particles and interactions: the term (^^) is needed for the renormalizability 
of the SM in the curved space-time. Since only one new parameter is introduced (the value 
of the non-minimal coupling ^, fixed by the amplitude of scalar perturbations at the COBE 
scale Q ) the Higgs inflation is a predictive theory, allowing to fix a number of parameters 
that can be found by the observations of the CMB. Moreover, the model does not have 
a problem of the graceful exit from the inflationary regime and leads to high values of 
the reheating temperature, Tr ~ 10^^ GeV |14]. The value of ^, necessary for successful 
inflation, is required to be rather large, ^ ~ 10^ — 10^. 

Due to the minimal character of the Higgs inflation it can be ruled out by the CMB 



observations and by particle physics experiments, especially those at the LHC. In refs. |15 



16] we conjectured that for the Higgs boson to play a role of the inflaton the SM has to 



be a consistent field theory all the way up to the Planck scale. This puts the mass of the 
SM Higgs to a specific interval. The mass is constrained from above by the requirement 
that no Landau pole appears in the scalar self-coupling of the Higgs boson below the 



Planck scale (see |T7| , |l8|, 19, ^] and references therein). The lower limit comes from the 
demand that the electroweak vacuum is stable for all scalar field values below Mp (see 
lH, m, m, m, m, m and refer ences therein) Our results were challenged in , claiming 
that radiative corrections to the tree Higgs potential are large in the inflationary region and 
that the spectral index Ug of scalar fluctuations is only consistent with WMAP observations 
for values of the Higgs mass niH ^ 230 GeV. The latter conclusions were shown to be 
modified drastically by the running of the coupling constants in [j^ and containing 
the one-loop and two-loop renormalization group (RG) improved analysis respectively. A 
similar remark was also made in |]28| . 

The aim of the present work is to upgrade the computation of [Q] to the two-loop level. 
This allows to check the stability of our conclusions against radiative corrections and to 
make a proper comparison with Q, where two-loop computation was presented. 

The paper is organised as follows. In section ^ we discuss the basic assumptions and 
elaborate on the possible choices of normalization point for radiative corrections. In section 
m we consider the action of the Standard Model in the inflationary region, in section ^ we 
compute the two-loop effective potential in the inflationary region. In section |5| we discuss 
the RG improvement of the effective potential. In section |6| we fix the procedure for 
computing the radiative corrections to the inflationary potential. In section we present 
the numerical results. In section |8| we discuss the difference between our analysis and that 
of ref. 1^. Section ^ is conclusions. We use the same notations as in Q. The present work 
can be considered as a continuation of this paper. In particular, we omit all the discussions 
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already made in [p. 



2. The basic assumptions 

It is a challenge to invent a self-consistent framework for computation of radiative cor- 
rections to inflation in any type of field theory. The reason is associated with non- 
renormalizable character of gravity. To clarify this point, let us consider the most popular 



line of reasoning which is based on the following logic (see, e.g. |2^, ^]), related to con- 
struction of effective field theories. 

Take some field theory coupled to gravity and consider cross-sections of different reac- 
tions (such as scalar-scalar or graviton-graviton scattering) computed in the lowest order 
of perturbation theory. Find the lowest energy A at which one of these cross-sections hits 
the unitarity bound. Call A the "ultraviolet cutoff". The fact that perturbation theory 
breaks down for energies above A may then signal that the theory under consideration is 
not a fundamental theory, but an effective one, valid only at momenta smaller than the 
ultraviolet cutoff A.^ If true, the initial Lagrangian must be modified by adding to it all 
sorts of higher-dimensional operators, suppressed by the scale A. Now, if the addition of 
these higher dimensional operators with coefficients of the order of one spoils the effect 
under consideration, it is said that the realisation of a particular phenomenon (inflation in 
our case) is "unnatural". At the same time, if one finds some symmetry principle which 
keeps the unwanted contributions of higher dimensional operators small or zero, then it is 
said that the corresponding effect is "natural".^ 

For the values of the non-minimal couplings ^ ^ 1 the "cutoff" scale, determined by 
the method described above, is of the order of:^ 

A. ~ ^ < Mp . (2.1) 

If the higher order operators, suppressed by A^, are added to the SM with coefficients of 
the order of one, they spoil the Higgs-inflation, and, according to the definition discussed 



above, make it "unnatural" |29, 3C]. In other words, the infiationary predictions depend 
in a sensitive way on the coefficients with which these operators appear. 

The point of view we take in this paper is very much different from the reasoning 
described above. Contrary to |2^, |3^, we will assume that the breaking of perturbation 
theory, constructed near the SM vacuum at energies greater than A, is not a signal of 
existence of new physics which should replace the SM, but rather an indication of a new 
physical phenomenon — strong coupling, which should be treated inside the SM by non- 
perturbative methods (such as resummation of different diagrams, lattice simulations, etc.). 

^This type of arguments were very successful for construction of renormalizable gauge theories with 
spontaneous symmetry breaking [^| in general and the SM in particular, starting from four-fermionic 
Fermi interaction. 

^Of course, the real physics question is not whether this or that theory is "natural" but whether it is 
realised in Nature. The small value of the scale of the electroweak symmetry breaking in comparison with 
the Planck scale is considered to be extremely "unnatural". The same is true for the cosmological constant. 
Nevertheless, these are the values we have to live with. 

^This estimate was done by S. Sibiryakov (private communication) and later in ji^, |3C| ]. 
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Put it in other words, we will assume that the SM is valid for all momenta smaller that 
the Planck mass Mp. Therefore, no higher dimensional operators, suppressed by will 
be added"^. 

Fortunately, the presence of the strong coupling at energies higher than in scat- 
tering of particles defined as excitations above the SM ground state (with small vacuum 
expectation value of the Higgs field) does not prevent the reliable computation of the in- 
flationary potential. As we will see (section in the inflationary region, for the Higgs 
field values /linf ~ Mp/y^, an adequate description of the system is given by the so-called 
chiral electroweak theory (SM with the frozen radial Higgs mode). This theory is 
not renormalizable and has an ultraviolet cutoff Ai^f ~ h^nfi defined with the use of the 
procedure, discussed above. Ai^f is associated with the scattering of the longitudinal inter- 
mediate vector bosons.^ At the same time, the typical energy saturating the loop diagrams 
appearing in the computation of the effective potential is of the order of the masses of the 
SM particles in the Higgs background (dimensional regularisation is assumed) and thus is 
smaller than the cutoff by the magnitude of gauge or Yukawa coupling constants. In other 
words, the contribution of the strongly interacting domain of energies to the low-energy 
observables such as the effective potential is suppressed. 

In spite of the fact that we do not include higher-order operators to the computation, 
the effective potential cannot be fixed unambiguously. In section 3 of [Q] we discussed two 
distinct possibilities. In the first one (referred as prescription I) the normalization point 
/X (coinciding, for example, with t'Hooft-Veltman parameter in dimensional regularisation) 
of the loop integrals is chosen to be proportional to the coefficient in front of the scalar 
curvature. In the Jordan frame this gives 

/x^ oc M|, + ^h^ , (2.2) 

and, consistently, in the Einstein frame 

^? (xMI . (2.3) 

With these choices the computations in Jordan and Einstein frames are equivalent both 
at classical and quantum levels. The prescription II corresponds to the Jordan frame nor- 
malization point given by (2.3). For consistency, one has to choose then the normalization 



point in the Einstein frame as 

(X Ml,/{M^ + ^h^) . (2.4) 

In short, the prescription I is "standard" (field-independent) in the Einstein frame, whereas 
the prescription II is "standard" in the Jordan frame. 

A possible physics behind the choice II is related to an idea that the Jordan frame is the 



one in which "distances are measured" [27| and that in this frame there exist a fundamental 



■^As has been shown in the higher order operators suppressed by Mp are harmless for inflation. 

^Not surprisingly, the energy at which the tree cross-sections start to exceed the unitarity bound, depends 
on the background field. For example, for the graviton-graviton scattering the corresponding cutoff is 
[Mff, rather than Mj,. 
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cut-off, independent on background scaiar fieid. A possible physics leading to the choice I 
is related to quantum scale invariance, discussed in |Q, though the prescription II can be 
easily realised in scale-invariant theories as well.^ 

To elucidate the difference between these two prescriptions let us consider the limit 
of large Higgs fields, ^/i^ S> Mp. In this case the dimensionful parameters -/VfJ, and the 
Higgs mass can be neglected, and the Lagrangian of the electroweak theory with grav- 
ity incorporated acquires a new symmetry (at the classical level) — the scale invariance, 
(/>(x) — > A"(/>(A3;), where ^ is a generic notation for any field with canonical mass dimension 
GeV^. Note that in the presence of gravity the scale invariance does not ensure conformal 
invariance, since C 7^ The fate of the dilatation symmetry at the quantum level depends 
on the way the infinities arising in loop computations are regularized and subtracted. 

The use of standard regularization and renormalization schemes (such as dimensional 
or Pauli-Villars) leads to the trace anomaly in energy-momentum tensor and thus to the 
breaking of dilatation symmetry (for a review see ||3^). The flat direction, existing for 
the Higgs field in Einstein frame, gets bended owing to radiative corrections. This is the 
prescription II. If gravity is taken away, it leads to a renormalizable field theory. 

A non-conventional approach is to use dimensional regularization (n = 4 — 2e) and keep 
the scale-invariance intact even if e / [^, |3^ (see also [^]) by replacing the t'Hooft- 
Veltman parameter // by a combination of dynamical fields with matching mass dimension. 
In this case the theory stays scale-invariant at the quantum level and the flat direction 
remains intact. This is the prescription I. If gravity is taken away, it leads to an effective 
field theory, valid for momenta transfers up to the effective Planck scale ~ [0^ ]. 

It is the trace (or conformal) anomaly in the total matter energy-momentum tensor 
which leads to the dependence of the inflationary predictions on the Higgs mass. The one- 
and and two-loop corrections to the effective potential, that we find in this paper in the 
inflationary region, account exactly for this anomaly, which is different in prescriptions I 
and II. 

It looks impossible to fix the unique prescription without knowing the physics at the 
Planck scales. In other words, the determination of inflationary cosmological parameters, 
such as the spectral index Ug through parameters of the SM, is subject to an intrinsic 
uncertainty, related to renormalization procedure. We will do computations with both 
prescriptions. Fortunately, the numerical difference happens to be small in two-loop ap- 
proximation, as it was so at one-loop [p. 

3. Action in the inflationary region 

After transformation to the Einstein frame the electroweak Lagrangian is essentially non- 
linear and is given by 

= \id^.xf - - -^H2 - L^/z + Ly- U{x) , (3.1) 

®In the language of the prescription I corresponds to the choice of the dimensional regularisation 
parameter /i^ oc ^xX^ +S,hh^ , while prescription II corresponds to the choice /i^ oc x^, where x is the dilaton 
field. 
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where 



Lw/z = -rtr:[V^] , 



Ly = iQL,RlPQL,R - {^QlUQr + ■■■ + h.c.) , 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



and D is the standard covariant derivative for the fermions. The field % is related to the 
module of the Higgs field h through the equation 



dx ^ 
dh 

with the confer mal factor Q.{h) given by 



^2 + Q(2i^2 



(3.6) 



(3.7) 



For h <^ Mp I the relation between % and h is 



1 + 6^2;^7M2 + 



arcsinh [\fQ^h/Mp) 



and for h S> Mp/^ we have 



Xc^\I^Mplog^}\h) . 



The parameter v in the inflationary case is 

v\h) 



(3.8) 



(3.9) 



(3.10) 



The dimensionless Nambu-Goldstone bosons n"' are parametrized in the non-linear form as 



U = exp [2i7r"r"] , = 



(3.11) 



In the Yukawa Lagrangian we only write explicitly the top quark contribution, and define 
the quark fields as 



Ql = {h, bL),U = -T2WT2 , Qr = (tp, 0) , 

where Ta are the Pauli matrices. 

The gauge fields and the fields strength are given by 



(3.12) 



B 



Y 



B^y — d^Bi, — d„Bn . 



(3.13) 
(3.14) 
(3.15) 
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The scalar potential is 



Uix) = ^ , (3.16) 



where we have neglected the Higgs mass, irrelevant for analysis of inflation. For h ^ Mp/^ 
we have 

f^W^l^(l-e"^""J • (3-17) 

As has been discussed in [Q], in the background of the "small" Higgs fields h < the 
action in the Einstein frame coincides with that of the canonical Standard Model. Thus, 
in this region the computation of the effective potential is straightforward (and in fact not 
necessary, as the inflation takes place at higher values of the scalar field). 

In the inflationary region the Higgs field is of the order 

> ^ , (3.18) 

corresponding to x ^ Mp. For ^ ^ 1 the Higgs field is essentially massless in the region 
(|3.18 ), and decouples from all the fields of the SM. The transition between these two regimes 



is governed by the effective couplings of the deviation 5x of the field % from its background 
value Xc to the fields of the SM. The canonical values of the Higgs field couplings to other 
fields are multiplied by the factor q given by 

.-^4 , ' (3.19) 

For small fields h <C Mp/^ the factor q is equal to the SM value q = 1 and for large 
fields Mp/^ h < Mpjyf^ it leads to the suppression of Higgs interactions by g = 
Mp/(^/6C/i) < 1. This effect can also be seen in the Jordan frame |2^, |2| , where it emerges 
due to the kinetic mixing between the Higgs field and gravitational degrees of freedom. 
The relation between the suppression factor s found in Q and the factor g is s = g^O^. 

In the limit /i — > co the action ( |3.1| ) is nothing but the chiral electroweak theory studied 
in detail in numerous papers (for a review see, e.g. ^^) plus a massless non- interacting 
field. It will be used for analysis of radiative corrections to the inflationary potential below. 



4. Two-loop effective potential 

The effective potential for the field x is gauge- invariant, since x is gauge- invariant itself. 
However, the computation is most convenient in Landau gauge, in which the Goldstone 
bosons (hidden in the matrix IX) and ghosts are massless, and a number of diagrams can 



be dropped due to the fact that the vector propagators are transverse |39]. The one-loop 
contribution to the potential is given by 

^m%r / mLr 5\ 3mi / mi 5\ 3mf / m? 3\ , ^ ^. 
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where we keep now the constant terms omitted in Eq. (9) of [Q], as they are needed for a 
consistent two- loop analysis. As in [||, the part associated with the (light in the inflationary 
region) field x is neglected. The background field masses of intermediate bosons [mw and 
mz)-, and of t-quark (rrit) are given in Eq. (8) of [Q]. 

To find the potential in two loops, one should first expand lA with respect to Goldstone 
fields tt". It is sufficient to keep first order terms in the Yukawa part ( |3.5| ) and second 
order terms in the gauge part (|3.1|). After this is done, the computation of the two- loop 
contribution is straightforward but rather involved. Fortunately, all necessary contributions 



can be extracted from |4C]. This paper found the following representation of the two-loop 
contribution: 

U2 = Vs + VsF + Vsv + Vfv + Vv , (4.2) 

where Vs represents scalar loops, Vsf — the diagrams involving scalars and fermions, Vsv — 
scalars and vectors, Vpv — fermions and vectors, and Vy — vectors and ghosts. To get the 
effective potential for our case one has to take away all diagrams involving the Higgs field, 
since the interactions of it with itself and all matter fields are suppressed by ^ in the 
inflationary region. In addition, all Goldstone masses have to be put to zero. As a result, 
the following substitutions have to be made: 



Vsf mi ( mflimt, 0,0) + ^J{mt,mt) 



Vsv 



[1 - 2sm^i 
2 cos2 e 



-A{0,0,mz) + A{0,0,mw) 



g sin 6rnzB{rnz, mw, 0) — e rnY^B{mw, 0,0) 



(4.3) 
(4.4) 

(4.5) 



whereas Vfv and Vv stay intact. Explicit expressions for different functions appearing 



here are quite long and we do not reproduce them, see |4C] for all definitions. 

The two-loop effective potential derived in this way contains an explicit dependence 
on the normalization point ^. It is, however, spurious as it must be compensated by 
the running of the coupling constants and the field x- If we knew the RG running of 
the parameters, the best choice of /i would be the background mass of a vector boson 
or of t-quark, which will minimize the logarithms in the higher-loop corrections. In the 
next section we will derive the necessary RG equations, valid for large field values in the 
inflationary region. 



5. Renormalization group equations 

The change of the scalar field from the border-line of applicability of the canonical Standard 
Model h ~ Mp/^ to the inflationary region /li^f ~ Mp/^/^ is relatively small. This means 
that the one-loop renormalization improvement of the two-loop effective potential will do 
already a very good job. The main objects which appear in the one-loop effective potential 
( |4.1| ) are the background masses of W , Z bosons and of t-quark. Therefore, we need to flnd 
the logarithmic running of these masses. It is sufficient to consider their renormalization in 
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the asymptotic region of large h, where the tree potential for the scalar field is flat. Here 
the action (3.1) is just the action of the chiral SM with v = Mpj^J^. 

To find the running of mz-, mw and mt it is sufficient to consider the WW, ZZ and tt 
two-point functions in the chiral SM. The explicit computation gives (in the MS subtraction 
scheme) : 

.2.. 9 2 _ (^9'^ 23^2 ^ 2 , %'^/ 



167r = QVt + -Y^j "^W > (5-1) 

16.V^™l = (^ + ^ + ^ - ^yi,f - - ^iy- + ^) ^ , (5.2) 

2 d 2 ( . .a 2\ 2 



167rV— m^^ = - + ml , (5.3) 

where nj = 3 is the number of fermionic generations. 

In fact, it is more convenient to rewrite these equations by introducing the running of 
the coupling constants g' , g, g^, yt, A, and ^. As the chiral SM is well studied in the liter- 
ature, we can extract the necessary results from the existing computations. Reference |41| 
gives the following equations for the running g' , g: 

^2,.d , _ (I 1 , 20n/\ ,3 



IQ^^n^g' =\ + L g'^ (5.4) 

.2„5_ ^43 , 1 4n;\ 3 



The number ^ accounts for the fact that the Higgs field (but not Nambu-Goldstone bosons) 
decouples from the fields of the SM. 

Reference also gives an equation for running of the parameter v"^ , without account- 
ing for the Yukawa couplings. It is not difficult to add the corresponding contribution, what 
leads to: 

16^V|;^' = + - QVt^ . (5.6) 

Since oc | in our case, this can be converted to equation for ^: 



d . /3 



IGvrV^e = - + - Qyt) e • (5.7) 

Clearly, the running of the strong coupling does not change in the chiral phase, so that 

167rV|^53 = -7gl . (5.8) 



Now, combination of Eq. ( p. 31) with (p. 61) gives an equation for the top Yukawa coupling 



2 9 / 17 ,0 3 



levr^^u— = [-r^g" " ^9^ " 893' + 3yt j yt ■ (5.9) 

The equation for A follows from equations for f ^ and from the fact that the one-loop 
effective potential U{h, A(/i)) -|- Ui{h,fi) in the asymptotic region of the Higgs fields must 
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be independent on /i: 

16vrV|;A = -6yt + ^ {2g' + (g" + g'f) + (-3^'^ - &g^ + I2yj) A . (5.10) 

This equation is quite peculiar as it does not contain the term and thus does not have 
the Landau-pole behaviour at one loop. It is instructive to provide also the equation for 
the combination A/^^, which is the only combination of constants present in the tree level 
potential ( p7|) 

This equation does not contain terms proportional to A at all, making the running less 
important for heavier Higgs masses. It is interesting to note that the right-hand-side of 
Eq. ( ^.11 ) is exactly the same for h < Mp/^, what can be derived from eqns. (10-14) of 
i- 

However, this system of RG equations is not complete. The reason is that the chiral 
SM is not renormalizable, so that in order to remove the divergences one has to add a 
certain set of counter-terms with the structures different from those already present in 
( |3.1| ). These operators were found and discussed in detail in a number of papers (see, e.g. 
a recent work ^^). Only few of these operators are needed for our purposes, namely those 
which contribute to the renormalization of the mass of Z boson. These are:^ 

-^aov^ {iT[TV^]f + ^QiB^,tr[TVFn , 

where T = 2L{T^U^ . The coefficients ao and ai obey the RG equations [41| 

2 d 3 ,2 2^1 

IGvr ^i— ao = -g , 167r fi—ai = - (5.12) 



and the mass of the Z-boson is 

= 7 ((5" + 9')(1 - 2ao) - 2g'^g^ai) v"" . (5.13) 



This expression, supplemented by the RG equations (^^, (^^, (|5l8|), (|5?7|), (|5^), ( |5.10| ), 



(5.12), allows for a one-loop improvement of the 2-loop effective potential in the inflationary 
region. The initial conditions for the running are to be fixed at the borderline between the 
SM phase and the SM chiral phase, corresponding to ^ ~ Mp/S^. 



6. The procedure for computations of inflationary parameters 

To determine the inflationary parameters, we generalize the procedure described in |Q] to 
two loop level. It consists of several steps described below. 

'^Note that the renormahzation of W and t masses does not require any new counter-terms, at least at 
the order in the coupUng constants we are working with. 

*A potentially important operator corresponding to ag (see for definition) does not contribute at 
the order in the coupling constants we are interested in. 
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To find the initial conditions for the RG running of the couphng constants at /i = 
we relate the values of the MS-scheme g, A and yt to the pole Higgs and top masses 
with the use of the formulas given in section 5 of ref. |42| and in the appendix of 



ref. 1 43]. The procedures described in these works give almost identical results (see 



appendix) . 

2. We replace the system of one-loop RG equations (10-14) of for all three gauge 
couplings, scalar self-coupling and top-quark Yukawa couplings by the corresponding 
two-loop equations taken from ||4^ and reproduced in convenient form in 1^^. Solv- 
ing these equations up to the scale Mp/^ gives us the values of parameters at the 
borderline between the SM and chiral SM. 

3. The values of the coupling constants found in the previous step are used as initial 
values for RG group running in the chiral phase of the EW theory. Since the param- 
eters uq and ai are absent in the SM, the initial conditions for them are fixed to be 

Q 

zero. 

4. The two-loop effective potential is computed in the inflationary region as U + U1 + U2, 
with all coupling constants taken at scale /x, and the Z boson mass given by ( 5.13| ). 



To optimize the amplitude of higher order corrections, the parameter /x is chosen as 
inQ. 

For the renormalization prescription I the value of /i is 

= K^^l{x) = -'^^ (1 - ^'^) ■ (6-1) 

Here k is some constant of order one.^'^ 

For the prescription II we take another value for /i: 

= mUxMxf = - 1) . (6.2) 

As the potential is fixed, the value of ^ is chosen in such a way that the correct 
WMAP5 normalization is reproduced. It is important to stress here that the wave- 
function normalization for the Higgs field in inflationary domain can be neglected, 
since the Higgs interactions with matter are suppressed by large value of non-minimal 
coupling ^. 
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The interpolation between the SM running of the coupling constants and the running in the chiral SM 



can be made smooth with the use of g-factor, defined in (3.1£), in analogy with ref. [g]. This leads, however, 
to negligible modifications of numerics, since q{x) defined in Eq. (3.19) is a rapidly changing function in 
the region h ~ Mp /^. A simple demonstration of this fact can be made by switching from the SM RG 
running to the chiral SM one at Mp/y^ instead of Mp/^. 



^°For technical reasons instead of expressions (6.1), ( |6.2| ) we used for numerics the same formulas, but 
with the constants yt, ^ taken at the point Mp/£^. 
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6. The spectral index ng and tensor-to scalar ratio r are computed. To be more specific, 
for the given potential the field x corresponding to the end of inflation is found, then 
the field corresponding to 59 e-foldings of inflation is found, then at this scale Ug and 
r are calculated using standard slow-roll formulas. 

The downgrade of this procedure to tree mapping at /i = mz, one-loop running from 
this point to Mp/^, and the use of one-loop effective potential in the inflationary region 
without RG improvement coincides with the computation carried out in our earlier work 
up to higher-order terms in coupling constants, due to Eq. (18) of Q.^^ The computation 



of 1 27] corresponds to the same procedure with the prescription II but without RG running 



from mz to Mp/^. In [27| the values of the coupling constants at /x = Mp/^ were taken 
to be the same as at /i = mz- The computation Q corresponds nominally (for a more 
detailed comparison see section ^) to one-loop mapping at /x = and the use of tree 
effective potential improved with the use of two-loop RG equations up to the inflation 
scale. 



7. Numerical results 

The results of the computation of the spectral index and tensor-to-scalar ratio using the 
procedure from the previous section with both prescriptions I and II are shown in Figs. |l|, 
g The running of the spectral index dns/dln k is always very small, of the order of 10 . 

Let us first comment on the results in the prescription I. The spectral index is nearly 
constant, exactly the same as for one- loop computation [Q]. The only difference is the shift 
of the minimal allowed Higgs boson masses to: 



mL 



mt - 171.2 as - 0.1176 

126.1 + — X 4.1 X 1.5 

2.1 0.002 



GeV , (7.1) 



which happens mainly due to the two-loop running up to the Mp/^ scale. One can see small 
wiggles at the lowest mn on the lines (they were more pronounced in the one-loop result). 
This corresponds to the situation when A is approaching zero in the inflationary region 
and the tree level part of the potential becomes comparable with the loop corrections. It 
is unclear whether this feature is a real physical effect or a result of approximation. 

At the high mu, contrary to the one- loop analysis, the solution to two- loop RG equa- 
tions for the scalar self-coupling does not show up a Landau-pole behaviour. Though A 
enters into strong-coupling region, it stays finite. Of course, perturbation theory cannot 
be trusted in this case. We believe that infiation is hardly possible in this situation. The 
lattice simulations of pure scalar theory [ 45 1 show that the corresponding theory has noth- 



ing in common with a continuum strongly coupled field theory. Requiring that one should 



^^The RG equations in the chiral phase of the EW theory found in the present work are different from the 
RG equations used in |l| in this domain. However, the one-loop running of the combination of couphngs 
A/^^, most relevant for determination of the inflationary potential, is exactly the same in both phases, see 
Eq. 
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Figure 1: Spectral index Us depending on the Higgs mass mu-, calculated with the RG enhanced 
effective potential. Nearly horizontal coloured stripes correspond to the normalization prescription 
I and different mj. Green, red, and blue stripes give the result with normalization prescription II 
for different mt and as = 0.1176, two white regions correspond to different as and mt = 171.2 GcV. 
The width of the stripes corresponds to changing the number of e-foldings between 58 and 60, or 
approximately one order of magnitude in reheating temperature. 
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Figure 2: Tensor-to-scalar ratio r depending on the Higgs mass mn, calculated with the RG 
enhanced effective potential. Nearly horizontal solid lines correspond to the normalization pre- 
scription I. Green, red, and blue dashed lines give the result with normalization prescription II for 
mt — 169.1, 171.2, 173.3 GeV. Dependence on the number of e-foldings is very small. 



have the weak coupling (somewhat arbitrary we set the requirement A < 6) in the region 
Mp/^ we obtain, in 2-loop approximation, 



rUy. 



193.9 + 



mt - 171.2 
2T 



X 0.6 



a, - 0.1176 



0.002 



X 0.1 



GeV 



(7.2) 



Note that the requirement to have weak coupling up to ;U = Mp is somewhat stronger, 
rriH < 173.5 GeV (for central values of mt and a^). Also, the requirement that A(/u) > 
for fi < Mp is just slightly stronger, than (|7.1|), giving the interval of the Higgs masses for 
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Figure 3: The running of tlie Higgs coupling constant A up to the high scale in the canonical SM. 
Solid line graphs correspond to mn — 174 GeV and itih — 126.3 GeV, leading to A = 6 at Mp, or 
A touching at a scale slightly below (these values represent the window in which the SM is valid 
all the way to the Planck scale, for central values of rrit and a^). Dashed graphs correspond to 
niH — 194 GeV and m^f = 126.1 GeV, leading to strong coupling or zero at Mp/^ (for ^ = 11600 
and 704, respectively). The lower part of the graph is scaled in the vertical direction for visibility. 




which the SM remains a good quantum field theory all the way up to the Planck scale, see 
Fig. I 

Figures § and ^ give for reference the values of ^ and A at the scale Mp/^. One can 
see, that the tree level COBE normalization requirement |14] holds very well for high Higgs 
masses, and becomes violated at lower ones, see Fig. 

For prescription II the behaviour is qualitatively similar to the one-loop calculation in 
m. Numerically, we have the overall shift in the Higgs mass region of 10 GeV, corresponding 
to the difference between one- and two-loop running of the coupling constants up to the 
inflationary region Mp/^. Also, the deviation from the asymptotic value of the spectral 
index Ug at small Higgs masses is smaller, than given by the one loop calculation. For low 
Higgs masses the spectral index drops, and goes out of the allowed region {rig > 0.93) for 



1 1 fcvmr 



mi -171.2 0^-0.1176 

126.7 + — X 4.5 X 1.7 

2.1 0.002 



GeV . (7.3) 



This is a bit stronger, than ([7.1|). 

It is interesting to note, that at large Higgs masses in any calculation the spectral 
index approaches the tree level value Ug = 0.967. The reason is that radiative corrections 
play a subdominant role if A is not very small, and the running of the coefficient A/.^^ in 
front of the tree level potential is less important for large A (see Eq. (|5.1l|) ). 
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Figure 4: ^ at the scale Mp/^ depend- 
ing on the Higgs mass mn for mt = 
169.1, 171.2, 173.3 GeV (from upper to lower 
graph). Solid lines correspond to prescrip- 
tion I, dashed — to prescription II. Chang- 
ing the e-foldings number and error in the 
WMAP normalization measurement intro- 
duce changes invisible on the graph. 
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Figure 5: A at the scale Mp/^ depend- 
ing on the Higgs mass mn for rrit = 
169.1, 171.2, 173.3 GeV (from upper to lower 
graph). Changing the e-foldings number 
and error in the WMAP normalization mea- 
surement introduce changes invisible on the 
graph. 



The comparison between different approx- 
imations for the effective potential can be 
found in Fig. 0. For this figure we always 
made two-loop SM running of the constants 
up to the Mp/^ scale, to evade the overall 
mass shift and allow for comparison. Then, 
we plotted one-loop effective potential ([4.1| ) 
and two-loop potential (^]^) without any fur- 
ther RG improvement; tree level potential 
with constants improved by chiral RG run- 
ning up to the scale ( |6.1| ) (or ( |6.2| ) for pre- 
scription II) , and RG enhanced tree level po- 
tential plus one- or two-loop potential. 

In prescription I all results are almost 
indistinguishable, because the normalization 
scale ( |6.1| ) changes very little during the in- 
flation. 
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Figure 6: The ratio ^/VA/46900 at the scale 
Mp/^. The ratio equal to one |lj] corresponds 
to the COBE normalization requirement for tree 
level potential. Solid lines are for prescription I, 
dashed — prescription II. Green, red and blue 
lines (from left to right) correspond to mt — 
169.1, 171.2, 173.3 GeV. 



In prescription II the results differ. We can see, that addition of the two-loop cor- 
rections changes the one-loop result rather strongly, because of accidental cancellation of 
the one-loop contributions from the gauge bosons and the top quark at one-loop. The 
RG enhanced result is rather close to the result obtained by using the two-loop effective 
potential. Furthermore, adding loop corrections to the RG enhanced tree-level result does 
not change the results much. 

Our results are stable against variation of the parameter k, defined in ( |6.lD , (|6.2|), 
demonstrating that the running of the coupling constants is indeed compensated by the 
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logs in the effective potential (as it should be). In Fig. |^ we show the change of the different 
approximations to effective potential due to the change of n. It is seen, that adding one- 
and two-loop logarithms to the RG-enhanced tree level potential removes the superficial 
dependence of the potential. 

The close similarity of the one-loop and two-loop results demonstrates that our con- 
clusions are stable against radiative corrections for both renormalization procedures. The 
main effect of accounting of two-loop running and one-loop matching is the shift of the 
admissible interval for Higgs boson masses. However, the spectral index Ug cannot be 
determined unambiguously in the SM inflation by giving the exact values of the SM pa- 
rameters (in particular the Higgs and top masses). Though remains close to its tree value 
in the wide range of Higgs masses, (small) deviations from it contain an extra uncertainty, 
coming from the choice of the normalization point (cf. prescriptions I and II discussed in 
section P). From theoretical point of view, this uncertainty can only be removed if the 
physics at Planck scales is fixed. A more exact cosmological measurement of ris, together 
with precise knowledge of the SM parameters can help to resolve this issue. 

It is interesting to note that the requirement of sucessful Higgs-inflation allows to 
put an upper bound on the mass of the top quark. If t-quark mass were larger than 
240-250 GeV, no choice of ^ and of the Higgs mass could lead to acceptable inflationary 
parameters. This number is smaller than the experimental lower constraint on the mass 
of the t'-quark of the fourth generation [46|, meaning that Higgs-inflation in the SM is not 
consistent with fourth fermionic family. 

We conclude this section with an estimate of the theoretical uncertainties related to 
higher order loop corrections^^. To test the accuracy of the pole matching procedure 
at low energy scale one can use somewhat different procedures to compute the values of 
the coupling constants at fx = mt- First, we simply used the formulas given in [|^, 43 1 
for fj, = mt. Then we used the same equations from |4^ at /i = mz, got the coupling 
constants at this point, and then considered their running up to /i = mt. The difference 
in the values gives an estimate of the uncertainties related to the two-loop terms in the 
matching procedure. For small Higgs masses in the region rriH — 130 GeV the variation 
in the scalar self-coupling is^'^ 5X{mt) / \{mt) ^ 0.018 and in the top Yukawa coupling 
^Utii^t) / UtiiT^t) — 0.0015 (uncertainties in other couplings are much smaller). This leads 
to the effective redefinition of the Higgs mass by 1.2 GeV and the top mass by 0.3 GeV. 
The change of the top mass by this amount leads to uncertainty in the lower limit on 
the Higgs mass equal to 0.6 GeV. The three-loop correction to the top mass gives 
5yt{fnt) / UtiTnt) — 0.0042 |^, the non-perturbative QCD effects in the top pole mass 
- MS mass matching are expected to be at the level of 5yt{mt)/yt{'mt) ^ 0.001 The 
four-loop as contribution to the top mass was guessed to be the same |5^. If we assume 
that these uncertainties are not correlated, and symmetric, we get a theoretical error in 
the determination of the critical Higgs mass, (Jmtheor — 2.2 GeV. This error is smaller than 
the present uncertainties related to the experimental errors in the top mass and in the 
strong coupling constant. It can be reduced further to 5mtheor — 0.4 GeV by the complete 
^^We thank Mikhail Kalmykov for discussion of this issue. 

^^We give two significant digits for uncertainties in A and yt in order to allow for a detailed comparison. 
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two- loop pole matching (not known at the moment), accounting for 3- loop as correction 
to yt, found in |47, |4^, and by the three- loop RG running up to the scale Mp/^ (note that 
the complete 3- loop /^-functions for the Standard Model are not known yet). 




U (5Mp ) /Mp 



130 140 150 160 170 180 190 
rtiH 

Figure 7: Spectral index depending on 
the Higgs mass mn calculated in different 
approximations. Nearly horizontal line rep- 
resents the result in the prescription I (dif- 
ferent approximations lead to indistinguish- 
able results). Everything else is the results 
in prescription II. Red dashed line (lower) is 
the one loop effective potential, green dashed 
line (upper) is the two loop effective poten- 
tial, close solid lines are RG enhanced re- 
sults without (blue) and with addition of 
the one or two loop effective potential (red 
and black, nearly indistinguishable). For all 
graphs mt = 171.2 GeV, a, = 0.1176. 
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Figure 8: Result of moving the normaliza- 



tion point (OA) in the potential. The black 
falling graph is the tree level contribution, 
red (lower) and green (upper) demonstrate 
removal of the dependence on the normal- 
ization point. Potential is calculated in pre- 
scription I for the typical inflationary field 
X = 5Mp and nin = 140 GeV. 



8. Comparison with refs. f^, ^ 

Our preprint [|[| appeared in hep-ph simultaneously with the paper devoted to the 
study of the same problem. The analysis of |^] was carried out with the use of two-loop 
RG equations for the running couplings, modified to account for non-minimal coupling of 
the Higgs field to Ricci scalar in the Jordan frame. In Q we used one-loop approximation 
and performed the computations of the effective potential in the Einstein frame. 

In very general terms, the findings of [|^, the present paper and of Q are quite similar. 
Namely, in |j^ it was shown that the Higgs-inflation in the SM, consistent with WMAP 
observations, is possible provided the Higgs mass mn lies in the interval mmin < "^H < 
iTimax, somewhat exceeding the region of the Higgs masses in which the SM is a viable 
effective field theory up to the Planck scale. In one-loop approximation the minimal and 
maximal Higgs masses were found to be 

mmin =[136.7 + {mt - 171.2) x 1.95] GeV , (8.1) 
mmax =[184.5 + {mt - 171.2) x 0.5] GeV . 
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The uncertainty resulting from variation of the strong couphng constant was not estimated 
in 0. 

With the use of two-loop analysis, the paper Q found that the successful SM inflation 
can take place if the Higgs mass exceeds the value 

m^in = [125.7 + I^LZHI X 3.8 - ~ "'^^^^ x 1.4] GeV ± 5 , (8.2) 
^ 2 0.0020 ^ ^ ' 

where 5 ~ 2 GeV indicates theoretical uncertainty from higher order corrections. 

At the quantitative level, the one-loop value of for mmin is some 11 GeV larger 
that that of [||. The result of our present paper for the minimal allowed Higgs mass is 
given in eqns. ( [7. 1[ ) , ( [7.3[ ) and almost coincides with the bound found in Q. Therefore, 
the shift of the critical Higgs mass is associated with the account for one-loop corrections 
at the low-energy matching (instead of tree procedure) and with two-loop (instead of the 
one-loop) running from the low scale to the high scale. The existence of this shift due to 
higher order effects was anticipated in |]l|. 

The behaviour of the spectral index as a function of the Higgs mass is a more subtle 



effect. The choice of the normalization point of Q coincides with that of [27| and corre- 
sponds to the prescription II of ||l| and of the present work. Let us compare our results for 
this case with those of [||. 

According to |^], if mn decreases, Us increases and is always larger than the tree value 
(see Fig. 1 of |||]). According to computations of ||l[ and of the present work, decreases 
and is always smaller than the tree value (Fig. 3 of and our Fig. ||). The similar type 
of discrepancy exists for the behaviour of the tensor-to-scalar ratio (cf. Fig. 5 of Q and 
our Fig. H). The two- loop computation carried out in the present work shows that the 
difference between the results of and Q is not related to the number of loops accounted 
for. Below we will try to elucidate possible origins of this discrepancy. Making a detailed 
comparison of computations in both articles is difficult as we used different frames for 
analysis and solved different equations. However, there are several essential points which 
are treated differently in these works even in one-loop approximation. 

(I) The first point is related to gauge invariance. Our method (and the one used in [l|) is 
explicitly gauge-invariant. The /3- functions for RG equations we used do not depend 
on the gauge choice, and the running of the couplings up to the point Mpj^ is gauge 
independent. The boundary conditions are gauge-invariant as well. The method of 
includes explicitly the gauge non- invariant object — the anomalous dimension of 
the scalar field 7. The equation (A. 5) of this work indicates that the Landau gauge 
was used (in the arbitrary a-gauge one would get 

7 = ^((2 + «)^^-3,?)+... (8.3) 

(here a is the coefficient in front of the longitudinal term of the gauge-field propagator, 
a = 1 and a = correspond to the Landau and Feynman gauges respectively). The 
gauge non-invariant parameter 7 is used then for RG running of all coupling constant 
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through equations hke dX/dt = + 7). Clearly, the definition of all couplings is 

then not gauge- invariant. It does not correspond to the standard gauge-invariant MS 



prescription. Therefore, the pole-MS matching described in [43| is not appropriate 
for fixing the initial condition for the gauge non-invariant couplings used in |0. The 
gauge variation of the scalar self-coupling A (most important for computation of 
inflationary spectral indexes) reads 

where 5a is the change of the gauge parameter. It does introduce the gauge depen- 
dence of the results at two-loop level and may result in the change of Us behaviour 
with the Higgs mass. 

(II) The second point is related to consistency of introducing of s-factor, which accounts 
for gravity-Higgs mixing in the Jordan frame. One can see that the argument leading 
to Eq. (4.6) of ||2| is only applicable to the physical Higgs field and does not work 
for Nambu-Goldstone bosons present in the Landau gauge, used in ||2|. Therefore, 
the kinetic term for these fields is not modified. In other words, the virtual Nambu- 
Goldstone bosons, contributing to diagrams for /3-functions must be kept. This fact 
was accounted for in our paper (section |5|) but not in as the comparison of the 
one-loop parts of relations (A.l) and (A. 2) and our equations in section ^ shows. 
One finds the differences in the running of the SU(2)xU(l) gauge couplings g and 
g' , of the scalar self-coupling A and of the top Yukawa coupling yt. For example, the 
/^-function for yt does not contain any contribution from yt in the inflationary region 
in whereas it is in fact present. We suspect that this may be one of the reasons 
for the opposite behaviour of the spectral index n<j found in ||2|, which underestimates 
this contribution. Moreover, the running of extra couplings oq and ai, necessary for 
description of the electroweak theory in the inflationary region, was not considered 
in§. 

(III) Yet another (numerically less important) difference is related to accuracy of compu- 
tation. Reference Q did not compute one- or two-loop effective potential but used 
the tree potential improved by the renormalization group. In this way the leading 
logs are accounted for, but the terms without logs are lost. In our present work all 
these terms up to the two-loop order are included. 

The second version of Q contains a "Note Added" in which a number of remarks about 
our work Q has been made. We comment on them below. For the reader convenience, we 
quote here the most relevant part of the "Note Added" of ref. Q : 

"In our analysis, we computed the full RG improved effective potential. We did this 
including (i) 2-loop beta functions, (ii) the effect of curvature in the RG equations (through 
the function s), (iii) wavefunction renormalization, and (iv) accurate specification of the 
initial conditions through proper pole matching. On the other hand, did not compute the 
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full effective potential or include any of the items (i)-(iv).^^ Instead ref. Q approximated 
the potential at leading log order with couplings evaluated at an inflationary scale after 
running them at 1-loop (this is one step beyond where couplings were not run)." 

We certainly agree with ||2[ that in |j^] we did not include (i) 2-loop (3 functions (we had 
one-loop running) and (iv) one-loop pole matching for initial conditions for renormalization 
group (we used the tree values). We said explicitly that our aim is the one-loop analysis 
of the problem, for which the tree pole matching is adequate. As for the point (ii), the 
one-loop renormalization group running of the combination A/^^, relevant for inflation, is 
exactly the same at small and at large Higgs field values, meaning that our analysis is 
perfectly consistent at one loop. Concerning the point (iii), our formalism does not require 
the computation of the wave-function renormalization. At the same time, the running of 
^ (which was not accounted for in Q) plays an important role. See also the discussion of 
wave- function renormalization and pole-matching procedure in (I) , and a comment on "full 
RG improved effective potential" in (III). 

9. Conclusions 

Inflation does not necessarily require the existence of the inflaton — the Higgs boson of 
the Standard Model can make the Universe flat, homogeneous and isotropic, and produce 
fluctuations necessary for structure formation. This happens in the SM with sufficiently 
large non-minimal coupling of the Higgs to the gravitational Ricci scalar. An important 
requirement is the validity of the SM all the way up to the Planck scale. In this case 
radiative corrections to the inflationary potential can be computed. They can be used to put 
constraints on the Higgs mass and to determine the connection between the cosmological 
parameters and properties of particles in the Standard Model. The window of allowed 
masses for the Higgs boson is niH S [126, 194] GeV (for central values of mt and and 
neglecting theoretical uncertaites). It somewhat exceeds the region where the canonical 
SM is a valid field theory up to the Planck scale, mn G [126.3, 174] GeV. This prediction 
can be testable at the LHC. The precision measurements of Us and r, together with exact 
knowledge of the SM parameters will allow to shed light on the origin of inflation and on 
Planck scale physics. 
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A. Initial conditions 

The MS values for the gauge coupling constants at // = mj were obtained from the /i = mz 

Though wavefunction renormahzation was not included in external leg corrections in the running 
of A were included. However these two effects roughly cancel against one another. 
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Figure 9: Correction to the Higgs self coupling Sh — Atrco — ^rrit) depending on the Higgs mass 
niH- 



present in |4g] by one- loop running for g' , g and two- loop running for gg. The resulting 
values for the gauge couplings at mt are g' = 0.3587, g = 0.6484, gg = 1.1619 + 4.5(qs — 
0.1176). Close values for electroweak constants were also obtained directly from the pole 



masses using formulas from [42|. 

The top quark Yukawa constant yt and Higgs self-interaction A at = can be 
obtained using the expressions from appendix of |43]. To summarize, the MS value for the 
top Yukawa is very well approximated as 



yt{mt) = ^^\ l + 5t) , 



rriH 



(A.l) 
130 GeV\ 



55000 GeV / 



(246.221 GeV)2. The initial data for A is \{mt) = mjj/2v'^ - 5h, where 



where v'^ = - 
6h is given in Fig. y. 

Using formulas from similar values are obtained. The main difference in 5t comes 
from the two- loop QCD contributions accounted for in |43] and not in The difference 
in 5h comes from neglecting the W and Z mass difference in the loop integrals in |42|, but 
accounted in |4^. Numerically, this just slightly changes the influence of the top mass on 
the overall results. 
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